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ABSTRACT 

The Magellanic Clouds are a local laboratory for understanding the evolution and properties of 
dwarf irregular galaxies. To reveal the extended structure and interaction history of the Magellanic 
Clouds we have undertaken a large-scale photometric and spectroscopic study of their stellar periphery 
(the MAgellanic Periphery Survey, MAPS). We present first MAPS results for the Small Magellanic 
Cloud (SMC): Washington M, T 2 + DD051 photometry reveals metal-poor red giant branch stars 
in the SMC that extend to large radii (~11 kpc), are distributed nearly azimuthally symmetrically 
(ellipticity=0.1), and are well-fitted by an exponential profile (out to i?~7.5°). An ~ 6 Gyr old, 
[Fe/H]«— 1.3 main-sequence turnoff is also evident to at least i?=7.3°, and as far as 8.4° in some 
directions. We find evidence for a "break" population beyond ~8 radial scalelengths having a very 
shallow radial density profile that could be either a bound stellar halo or a population of extratidal 
stars. The distribution of the intermediate stellar component (3 < R < 7.5°) contrasts with that of 
the inner stellar component (R < 3°), which is both more elliptical (e w 0.3) and offset from the center 
of the intermediate component by 0.59°, although both components share a similar radial exponential 
scale length. This offset is likely due to a perspective effect because stars on the eastern side of the 
SMC are closer on average than stars on the western side. This mapping of its outer stellar structures 
indicates that the SMC is more complex than previously thought. 

Subject headings: Galaxies: interactions — Local Group — Magellanic Clouds — Galaxies: dwarf — 
Galaxies: individual (SMC) — Galaxies: photometry 



1. INTRODUCTION 

The Magellanic Clouds (MCs) have had a complex his- 
tory that includes a recent close encounter with each 
other (about ^200 Myr ago), a presently ongoing pe- 
riod of intense star formation, and immense loss of gas 
to the various components of the Hi Magellanic System 
(the Magellanic Bridge, Magellanic Stream and Lead- 
ing Arm). Systematic study of the Small Magellanic 
Cloud (SMC) stellar morphology started with the pi- 
oneering work by de Vaucoulcurs (1955), who studied 
mostly the bright and young structures and noted the 
SMC's very irregular shape. Photographic plate pho- 
tometry of red horizontal branch stars (RHB) helped 
show that the older SMC population extends to a ra- 
dius of ^5°, is fairly azimuthally symmetric, but has a 
larger line-of-sight depth on the northeastern side than 
for the rest of the SMC likely due to the recent interaction 
with the Large Magellanic Cloud (LMC; Hatzidimitriou 
& Hawkins 1989; Gardiner & Hawkins 1991; Gardiner & 
Hatzidimitriou 1992). Additionally, Gardiner & Hawkins 
(1991) found that the SMC line-of-sight depth on the 
western side increases with radius, a finding more consis- 
tent with a spherically symmetric radial exponential law 
than with a flattened system (i.e., a disk). Stellar radial 
velocity studies (Kunkel ct al. 2000; Harris & Zaritsky 
2006) show no signs of rotation but instead that the stars 
arc pressure supported and very likely in a spherical or 
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ellipsoidal structure. In contrast, rotation is observed in 
the Hi component of the SMC for R <3°(Stanimirovic 
et al. 2004), which reinforces that the structure of the 
SMC is complex and multi-faceted. In deep CCD pho- 
tometry of three fields to the south of the SMC, Noel & 
Gallart (2007) found old SMC stars out to 5.8°, whereas 
De Propris et al. (2010) find very few SMC stars at a 
radius of ^5° toward the east and suggest the edge of 
the SMC is at a radius of 6 kpc. All of the above results 
for the SMC are generally in agreement with the follow- 
ing known properties of dlrr galaxies: the old and red 
stellar components (as traced by evolved stars) arc (1) 
spatially smoother and more extended than their younger 
and bluer counterparts, and (2) generally follow a radial 
exponential profile (e.g., Mateo 1998). However, many 
questions remain about the structure of the SMC pe- 
riphery beyond R^5°: Where is the edge of the SMC? 
Does the SMC have a stellar halo or stellar tidal tails? 
How has the interaction of the MCs with each other (and 
the MW) disturbed their stellar structures? 

We have undertaken the MAgellanic Periphery Survey 
(MAPS), a large photometric and spectroscopic survey 
of the stellar periphery of the MCs designed to elucidate 
their extended structure and kinematics and piece to- 
gether their past interactions as recorded in their outer 
morphologies. In this Letter, we present our first pho- 
tometric results for the SMC (§2). Using Washington 
M, T 2 + DD051 photometry we detect metal-poor SMC 
red giant branch stars to a radius of 10.6° (11.1 kpc) 
and at all of the many position angles explored around 
the SMC (§3). These stars are distributed nearly az- 
imuthally symmetrically and are well-fitted by a radial 
exponential profile to R w 7.5°. We find evidence for a 
"break" population beyond R <~ 8° with a shallower ra- 
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Fig. 1.— MAPS fields of the SMC periphery. Shallow fields are in 
red (M<22) while deeper fields are in purple (M<24). The central 
image shows the RGB starcounts from OGLE-III and MCPS. Color 
magnitude diagrams and color-color diagrams for the circled fields 
are shown in Fig. 2. 

dial density profile that could be a bound stellar halo or 
a population of extratidal stars ($ 4). The fitted center 
of the intermediate stellar distribution (3< R <7.5°) is 
offset by 0.59° from that of the inner stellar distribution, 
which is also much more elliptical This lack of centration 
is likely due to a perspective effect (§4). 

2. OBSERVATIONS AND DATA REDUCTION 

The stellar periphery of the SMC was imaged as part of 
the MAPS survey using the Washington M, T 2 + DD051 
photometric system (Majewski et al. 2000a), which sub- 
stantially helps discriminate foreground MW dwarf stars 
from the Magellanic red giant branch (RGB) stars. Ob- 
servations were obtained with the MOSAIC II camera 
(36'x36') on the CTIO-4m Blanco telescope on UT 2006 
February 28-March 4 and UT 2007 August 16-18 with 
the majority of the SMC data coming from the second 
observing run. In all fields net integrations of (40, 60, 
280)s in (M, T 2 , DD051) were acquired, and in a sub- 
set of 33 fields (most arranged in four radial "spokes" 
around the SMC) additional 300s observations were also 
obtained in M and T 2 (see Fig. 1). Observations were 
also obtained for Geisler (1990) standard star fields to 
enable photometric calibration. 

The SuperMACHO pipeline (Rest et al. 2005; Mik- 
naitis et al. 2007) was used to obtain flattened images 
and the PHOTRED pipeline (D. Nidever et al. 2011, in 
preparation) was used for the rest of the photometric 
reduction. PHOTRED makes use of a combination of 
the stand-alone DAOPHOT/ALLSTAR (Stetson 1987), 
ALLFRAME (Stetson 1994) and Sextractor (Bertin & 
Arnouts 1996) packages to make photometric measure- 
ments and create final, calibrated, dereddened photom- 
etry. Figure 1 shows our MAPS SMC fields and the 
density of SMC RGB stars in the central region (R<2°; 
color image) as selected from the OGLE-III (Udalski et 



al. 2008) and MCPS (Zaritsky et al. 2002) SMC pho- 
tometric catalogs. Examples of the final MAPS SMC 
photometry can be seen in Figure 2. 

3. THE RADIAL DENSITY PROFILE 

Figure 2 shows dereddened (M—T 2 , M) color magni- 
tude diagrams (CMDs) and (M-T 2 , M-DD051) color- 
color diagrams (2CDs) for three fields in the northern 
part of the SMC. These fields are typical for the north 
and northeastern SMC, which have higher densities than 
in other directions as described below. The SMC main- 
sequence (MS), RGB, and red clump (RC) are visible 
in the three fields extending out to i?=8.4°. The right- 
most panel is simulated photometry for the i?=8.4° field 
using the Besancon Galactic model (Robin et al. 2003) 
to ascertain the approximate distribution of MW fore- 
ground stars. Padova isochrones (Girardi et al. 2002) 
with [Fe/H]=-0.6/0.7 Gyr, [Fe/H]=-1.0/2 Gyr, and 
[Fe/H]=— 1.3/6 Gyr, shifted to a distance of 63 kpc, are 
also shown (red, gold and blue respectively) as fiducials. 
The dominant population in the left three panels is ^6 
Gyr old but a sizable younger population (~2 Gyr) is 
also visible to at least i?=6.2° as well as a weaker ^0.7 
Gyr population. Giant and dwarf stars separate well in 
the 2CD due to the gravity sensitivity of the DD051 fil- 
ter. Our giant selection criteria are shown by the dashed 
line and all stars selected as giants (with M o <21.0) are 
plotted as filled red circles. 

We derive starcounts in our fields by using a 2CD giant 
selection (dashed line in Fig. 2 lower panels) and an SMC 
RGB-CMD selection (dashed line in Fig. 2 upper panels). 
Our starcounts can be contaminated by MW halo giants 
and metal-poor MW dwarfs that fall in the 2CD giant 
selection. To estimate and remove this contamination we 
use the same 2CD selection criteria and the same CMD 
selection criteria but offset by —3 magnitudes in order 
to not overlap the SMC stars (following Majewski et al. 
2000b; see, e.g., their Fig. 12). This shifted CMD se- 
lection extends beyond our bright limit; to account for 
missing giants we use a related dataset (Majewski 2004) 
with the same photometric filters and reduced with the 
same software but that extends to brighter magnitudes 
in fields far enough from the MCs not to be contami- 
nated by them; from these we calculate that we are on 
average missing 22% of the foreground giants due to our 
brightness limit 4 . This constant fraction is used to cor- 
rect the foreground values for all fields. The selection 
at these brighter magnitudes should have contamination 
levels similar to those in our actual SMC selection if the 
Galactic halo field giants follow a ~i? -3 law, which pro- 
duces a uniform distribution in magnitude (Majewski et 
al. 2000b). This method has been used (and more ex- 
tensively explained) in several related studies (Majewski 
et al. 2000b; Palma et al. 2003; Westfall et al. 2006) and 
shown to be reliable (Majewski et al. 2005). Finally, we 
perform a 2D spatial linear fit to the foreground contam- 
ination levels of all fields to obtain a smooth distribution. 

The starcounts as a function of radial distance from 
the SMC center (taken as a=00 h 52 m 44 s , c5=-72°49.7'; 
Mateo 1998) can be seen in Figure 3a. The median fore- 
ground level is shown. Fields with R > 7.5° towards the 

4 This fraction is small because in general the RGB-tip is less 
well populated. 




Fig. 2.— Photometric data for three "deep" SMC MAPS fields at _R=4.0, 6.2 and 8.4° (circled in Fig. 1). The top panels show 
Schlegel et al. (1998) dereddened (M-T 2 ,M) CMDs with some fiducial Padova isochrones shifted to 63 kpc ([Fe/H]= -0.60/0.7 Gyr in 
red; [Fc/H]= —1.0/2 Gyr in gold; and [Fe/H]= —1.3/6 Gyr in blue). The dashed line shows the region used to select SMC RGB stars. 
The bottom panels show dereddened (M-T2, M-DD051) 2CDs for Mo <21.0 with the giant selection (dashed line). Giants are plotted 
with filled red circles in both top and bottom panels. SMC giants and MS stars are visible in all three fields. The right-most panel shows 
simulated photometry for the R=8A° field using the Besancon Galactic model to ascertain the approximate distribution of MW foreground 
stars, which are roughly bounded (on the blue edge) by the green line (also shown in the i?=8.4° panel). 



LMC (four fields) were removed because these are heav- 
ily "contaminated" by LMC RGB stars. We detect SMC 
stars above the foreground in all fields to i?=10.6° (11.1 
kpc), almost twice as far as the previous most distant 
photometric detection of SMC stars (R ~ 5.8°, Noel & 
Gallart 2007). 

The SMC starcount density follows a general radial 
exponential, but the distribution is not well-fitted by a 
single exponential for all radii and position angles, as 
can be seen by the large scatter about the exponen- 
tial fit (dashed line, fit to R < 8° fields). The density 
along each radial spoke follows roughly a radial expo- 
nential profile with nearly identical radial scale lengths 
(h r ss 0.8°) but with different amplitudes - northeast- 
ern fields having higher amplitudes than southwestern 
fields. A simple explanation of this pattern is that the 
center we have used - the center of the inner stellar dis- 
tribution - is not appropriate for the outer SMC stel- 
lar distribution. Figure 3b shows an elliptical exponen- 
tial fit to the same fields that also allowed the den- 
sity center to vary. This provides a much better solu- 
tion than the previous fit, with best-fitting parameters: 
accntcr=01 h 00 m 31 s , 5 C ontcr=-72°43'll"±4' (0.59° from 
the Mateo 1998 center), ^=1.0°, ellipticity=0.10, and 
line-of-nodes=154.4° (E of N). An illustration of the fit- 



ted distribution is shown in Figure 4. The shape of the 
outer contours (R>7.5°) is not well constrained by our 
data since we only sample a few position angles. For Fig- 
ure 4 the ellipticity and line-of-nodes of the intermediate 
distribution were simply extended to the outer distribu- 
tion. 

The distribution of RGB stars is quite different for the 
inner (i?<3°) and intermediate (3<i?<7.5°) radial re- 
gions. The intermediate component is much more az- 
imuthally symmetric (e=0.1) than the inner component, 
which is quite elliptical (e=0.3; Harris & Zaritsky 2006), 
and the major axes of the two components are nearly 
perpendicular (difference of ~105°; inner line-of-nodes 
50°; Harris & Zaritsky 2006). The center of the interme- 
diate component is offset by 0.59° to the northeast of the 
center of the inner component (red and black crosses re- 
spectively in Fig. 4) - along the inner's major axis - and 
is closer to the Hi dynamical center (separation=0.38°; 
Stanimirovic et al. 2004). 

Finally, Figure 3b, and the sudden increase in contour 
spacing seen in the outer parts of Figure 4, both show 
that there is a striking departure from the exponential 
fit for fields at large radii (R>7.5°). This "break" popu- 
lation extends from i?«7.5° to our fields at i?=10.6° and 
has a much shallower radial density profile (/i r «7°). Be- 
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Fig. 3. — The SMC radial density profile, (a) Giant stareounts from our MAPS fields as a function of radial distance from the Mateo 
(1998) SMC center. Color indicates the position angle (PA) of the fields (East of North; PA=116° is towards the LMC) and lines connect 
the fields of the four "radial spokes" (for the inner three fields). The best-fit exponential profile (to the R <7° fields) with h T =0.8° is shown 
by the dashed line. The dotted line indicates the median foreground level. Large offsets in the stareounts are apparent and indicate that 
the center used is not appropriate for the SMC periphery. (6) Giant stareounts as a function of major axis radius for our best-fit elliptical 
exponential model. This profile is a much better fit to all of the data than the model in (a). A "break" population at R>7.5° is apparent. 
The estimated foreground contamination levels for each field are shown as open squares. 



cause the detected break population density is only a 
few times the foreground density, more data are needed 
to determine its density at these large radii securely. 

4. DISCUSSION 

It is common for dlrr galaxies to have patchy distri- 
butions of young, blue stars but smoother and more ex- 
tended distributions of older, redder stars that follow ex- 
ponential (or King) surface density profiles (Mateo 1998). 
The fairly azimuthally symmetric elliptical exponential 
profile of the older SMC stars is consistent with this 
general trend. However, this symmetrical distribution 
(e~0.1) of the stellar periphery is quite surprising in the 
specific case of the SMC considering (1) the SMC's re- 
cent close encounter with the LMC ^200 Myr ago, which 
stripped large amounts of Hi from the SMC to create the 
Magellanic Bridge (e.g., Muller & Bekki 2007), (2) that 
the MCs have just recently passed periGalacticon (e.g., 
Besla et al. 2007) and thus were subjected to stronger 
tidal forces from the MW, and (3) that the LMC's outer 
stellar contours are more elliptical (e«0.2; van der Marel 
2001) than the SMC's even though the LMC is lOx more 
massive. How is this possible? 

On closer inspection, we find that the SMC is not as 
"undisturbed" as might be suggested from the fairly sym- 
metrical density profile of its intermediate component. 
Hatzidimitriou et al. (1989) and Gardiner & Hawkins 
(1991) used red clump stars to show that the line-of-sight 
depth varies significantly across the face of the SMC. It is 
much thicker ( "puffed- up" ) in the northeast and extends 
to smaller distances there (from the Sun) compared to 
the southwest. Therefore, the three-dimensional struc- 
ture of the SMC is somewhat irregular and disturbed, 
likely due to the recent close interaction of the SMC with 
the LMC. The fact that the projected density profile is 
still so symmetric might be a coincidence owing to the 



particular orientation of the MC system relative to our 
viewing angle. 

The difference in the line-of-sight depths between the 
eastern and western sides of the SMC suggests a possible 
explanation for the lack of centration between the inner 
and outer distributions. Because the stars on the east- 
ern side of the SMC are on average closer than the stars 
on the western side, the SMC periphery will exhibit a 
perspective effect. The near-side of the SMC will appear 
larger than the far-side and the density contours will be 
"stretched" on the near-side (east) and "bunched-up" on 
the far-side (west). The net result is that the center of 
outer contours of a fit to the on-sky densities will be 
systematically shifted to the near-side, or the east, as 
observed. A perspective effect has previously been ob- 
served to affect the outer contours of the LMC due to 
the inclination of its stellar disk (van der Marel 2001). 
However, it is important to note that the lack of any stel- 
lar rotation seen in radial velocities (Kunkel et al. 2000; 
Harris & Zaritsky 2006) or proper motions (Piatek et al. 
2008) indicates that the SMC does not have a stellar disk 
(like the LMC), but is pressure supported and likely has 
a more spheroidal or ellipsoidal shape. The difference in 
line-of-sight distances across the SMC might be a result 
of tidal distortions of the SMC shape due to the recent 
encounter of the MCs ^200 Myr ago. 

What is the nature of the break population that ap- 
parently dominates over 7.5<i?<10.6° and is fairly az- 
imuthally symmetric for at least 270°? The two main 
possibilities are (1) a tidal tail or debris, or (2) a "classi- 
cal" bound stellar halo as exists around the MW and M31 
(e.g., Guhathakurta et al. 2005). The MCs had a recent 
close encounter (^200 Myr ago), which no doubt had 
a profound impact on the SMC. The LMC very likely 
stripped Hi gas from the SMC to form the Magellanic 
Bridge, and possibly also puffed-up the stellar distribu- 
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Fig. 4. — Stellar surface density of the SMC. The central col- 
ored image shows the RGB starcounts from OGLE-III and MCPS. 
Contours show our best-fit model to the SMC RGB starcounts in 
MAPS fields and are drawn at even intervals in log density (start- 
ing at 33 giants deg - 2 and increasing with an interval of 0.08 in 
log density). Squares show MAPS fields constraining the model 
fit. The optical center (black cross), HI dynamical center (Stan- 
imirovic et al. 2004, blue cross), and model center (red cross) are 
indicated. 

tion on the northeastern side of the SMC. It's possible 
that this encounter also created a stellar tidal stream. 
Extratidal stars energy sort into the well-known double- 
tidal tail shape only outside a few tidal radii and at 
smaller radii can appear fairly azimuthally symmetric 
(Munoz et al. 2008). Since the tidal radius of the SMC 
is 4-9 kpc (Stanimirovic et al. 2004) the break popula- 
tion is within ^2 tidal radii; it is, therefore, difficult to 
distinguish a bound classical halo from extratidal stars 
in this regime. An anlysis of the kinematics and velocity 
dispersion of these stars as well as further photometric 
mapping at more position angles and larger radii might 
help elucidate their true nature (i.e., bound versus un- 
bound). Either way, these stars can be thought of as a 
newfound "halo" component of the SMC. 

De Propris et al. (2010) claim that the edge of the 
SMC is at a radius of ~6 kpc on the eastern side because 
they detected only five spectroscopically-confirmed SMC 
giants in an eastern field at i?^5 kpc. We have a field 
~1.6° away (but at a similar radius) from theirs and 
find 57 giant candidates with our Washington+DD051 
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selection method (Af<21.0). When we also use the De 
Propis et al. 2MASS selection criteria (13 < K s < 14, 

0. 5 < J—K s < 1.5) this reduces the number of giants to 
five in agreement with the De Propris et al. result. The 
difference in detected giants between the two techniques 
is due to the shallowness of the De Propis et al. selection, 
which is limited to Af~T7 at the color of the SMC RGB 
(M-T 2 ~1.7). As can be seen in the i?=6.2° panel of Fig- 
ure 2 this selection only samples the very tip of the SMC 
RGB whereas our selection goes four magnitudes fainter 
and is more sensitive to lower SMC densities. Therefore, 
we find the De Propris et al. conclusion that the SMC 
edge toward the east is at i?~6 kpc to be premature; our 
deeper data show that the SMC extends to much larger 
radii (at least i?~9 kpc and likely to i?~ll kpc, Fig. 3). 

5. SUMMARY 

In this paper we have discovered SMC stars to a radius 
nearly twice as large as the previous most-distant SMC 
detections. We observe two outer SMC stellar structures: 

1. An intermediate component of older stars dominates 
over 3<i?<7.5°, follows a slightly elliptical exponen- 
tial profile (/i r =1.0°), is nearly azimuthally symmetric 
(e=0.1), and has a center that is offset from the center 
of the inner stellar distribution by 0.59° (to the north- 
east) likely due to perspective effects. The structure of 
this component is probably spheroidal or ellipsoidal. 

2. An outer component dominating 7.5<i?<10.6° that is 
fairly azimuthally symmetric over at least 270° and 
has a shallow radial scale length (ft r «7°). The com- 
ponent could be a bound stellar halo or a population 
of extratidal stars. Analysis of kinematics and veloc- 
ity dispersions of these stars and further photometric 
mapping are needed to better reveal the nature of this 
component. 

The SMC, therefore, appears to be much more complex 
than previously thought and composed of several struc- 
tural components, like larger galaxies. 

D.L.N, was supported by an Sloan Digital Sky Survey- 
Ill APOGEE software postdoc and the NSF grant AST- 
0807945. We acknowledge funding from NSF grants 
AST-0307851 and AST-0807945, and NASA/JPL con- 
tract 1228235. We thank the OGLE and MCPS projects 
for making their SMC photometric databases available 
to the public, and the anonymous referee for useful com- 
ments that improved the manuscript. 

Facilities: CTIO (MOSAIC II). 



Guhathakurta, P., Ostheimer, J. C, Gilbert, K. M., Rich, R. M., 
Majewski, S. R., Kalirai, J. S., Reitzel, D. B., & Patterson, R. J. 
2005, arXiv:astro-ph/0502366 

Harris, J., & Zaritsky, D. 2006, AJ, 131, 2514 

Hatzidimitriou, D., Hawkins, M. R. S., & Gyldenkerne, K. 1989, 

MNRAS, 241, 645 
Hatzidimitriou, D., & Hawkins, M. R. S. 1989, MNRAS, 241, 667 
Kunkel, W. E., Demers, S., & Irwin, M. J. 2000, AJ, 119, 2789 
Majewski, S. R., Ostheimer, J. C, Kunkel, W. E., & Patterson, 

R. J. 2000, AJ, 120, 2550 
Majewski, S. R., Ostheimer, J. C, Patterson, R. J., Kunkel, W. E., 

Johnston, K. V., & Geisler, D. 2000, AJ, 119, 760 



6 



NIDEVER ET AL. 



TABLE 1 
MAPS SMC Densities 



Field Name" 


RA 


DEC 


-Rsmc 


RSMC, ellipse 


PAsmc 


-^giants 


Nback 


PSMC 


PSMC Error 




(J2000) 


(J2000) 


(deg) 


(deg) 


(deg) 






stars deg -2 


stars deg -2 


190L225a 


01:09:37.5 


-75:05:20.2 


2.54 


2.47 


154.71 


2581 


9 


7,142.6 


141.4 


190L225b 


01:09:37.7 


-74:32:07.8 


2.08 


1.93 


147.21 


4588 


9 


12,718.0 


188.3 


40S026 


01:12:44.0 


-69:10:20.5 


4.00 


3.90 


26.38 


549 


8 


1,502.5 


65.6 


40S071 


01:39:58.5 


-71:10:10.1 


4.00 


3.82 


71.25 


664 


8 


1,819.6 


72.1 


40S116 


01:46:01.8 


-74:11:49.4 


4.01 


3.68 


116.21 


862 


9 


2,367.8 


82.0 


40S161 


01:15:02.7 
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18.6 
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17.0 
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10.03 


251.21 
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14.6 
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9.87 
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9.49 


10.70 


206.13 


32 


12 


54.8 


18.5 
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a Ficld names are of format "radius-L/S-PA" where radius is 10 X degrees, L=LMC center, S=SMC center, PA=position 
angle with respect to L/S. 
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